################ 5.1 Editando alguns parmetros grficos

setwd("C:/Users/Renan Parmigiani/Documents/MESTRADO/Disciplina R/Exercicio 5")
esaligna <- read.csv("esaligna.csv", header = TRUE, sep = ",", as.is = TRUE)
head(esaligna)
str(esaligna)
plot(esaligna$dap~esaligna$ht) ##Vendo o grfico sem os ARGs
par (cex = 1.4) 
plot(esaligna$dap~esaligna$ht, xlab = "Altura (m)", ylab = "DAP (cm)", main = "Relao entre DAP e altura", tcl = 0.5, bty = "l") 
dev.off() ## Retornando os valores de par para o default

############### 5.2 Dois grficos juntos

esaligna
esaligna$talhao <- factor(esaligna$talhao) 
str(esaligna$talhao)
class (esaligna$talhao)


#Agregando valores por talho
valores <- aggregate(esaligna$dap, by = list (esaligna$talhao), mean)
sd <- aggregate(esaligna$dap, by = list (esaligna$talhao), sd)
valores$sd <- sd$x
valores
colnames(valores) <- c("talhao", "media", "desvio.p")
valores$cima <- valores$media+valores$desvio.p
valores$baixo <- valores$media-valores$desvio.p


par(mfrow = c(1,2))
boxplot(esaligna$dap~esaligna$talhao, main = "a", xlab = "Talho", ylab = "DAP", ylim= c(0,25))
plot(valores$media~valores$talhao, main = "b", xlab = "Talho", ylab = "DAP", ylim= c(0,25))
segments(x0 = (1:6), y0 = valores$baixo, y1 = valores$cima)
dev.off()
  

###################### 5.3 Adivinhando o cdigo


dados <- read.csv("arquivo.csv", header = TRUE, as.is = TRUE, sep = ",")
dados
str(dados)
head(dados)


summary(dados) # Buscando valores mnimos  emximos para comparar com os grficos apresentados
dev.off()
par(mfrow= c(1,2), tcl = 0.4)
par(mar=c(5,5,5,1))
par(pch = 17, bty = "l")
plot(dados$y1~dados$x1, xlab = "Log(Patch size)(ha)", ylab = "", xlim = c(0.5,2.2), ylim = c(0,3))
mtext(side = 2, "Euclidean distances", cex = 1.3, padj = -2.5)
text(x = 2.1, y=3, "a", cex = 1.5)
linha<-lm(dados$y1~dados$x1)
segments(x0 = 0.5, y0 = 2.25, x1=2.1, y1=1, lwd = 2 )

par (bty = "l", tcl = 0.4)
par(mar=c(5,1,5,3))
boxplot(dados$x2~dados$y2, ylim = c(0,3), outline = FALSE, names = FALSE)
mtext(c("Small", "Medium\n Edge", "Medium\n Inferior", "Large\n Edge", "Large\n Inferior", "Control"),side= 1, cex=0.8, at = c(1:6), padj = 1)
text(x = c(1:5), y=3, c("*","*","**","*","***"), cex = 2)
text(x = 6.6, y = 3, "b",cex = 1.5)

